Cooling a quantum circuit via coupling to a multiqubit system 
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The cooling effects of a quantum LC circuit coupled inductively with an ensemble of artificial 
qubits are investigated. The particles may decay independently or collectively through their in- 
teraction with the environmental vacuum electromagnetic field reservoir. For appropriate bath 
temperatures and the resonator's quality factors, we demonstrate an effective cooling well below 
the thermal background. In particular, we found that for larger samples the cooling efficiency is 
better for independent qubits. However, the cooling process can be faster for collectively interacting 
particles. 
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I. INTRODUCTION 

The ability to cool interacting quantum systems be- 
low the values imposed by the thermal fluctuations of 
the environmental reservoir of each subsystem is actu- 
ally of great interest For instance, a laser cooling 
technique for trapped particles exploiting quantum inter- 
ference, or electromagnetically induced transparency, in a 
three-level atom was presented in Q ■ There, by appropri- 
ately designing the absorption profile with a strong cou- 
pling laser, the cooling transitions induced by a cooling 
laser are enhanced while heating by resonant absorption 
is strongly suppressed. The experimental demonstration 
of ground state laser cooling with electromagnetically in- 
duced transparency was performed in Q. Recently, this 
idea was adopted to cool a nanomechanical resonator 
0| • Further, the collective-emission- induced cooling of 
atoms in an optical cavity was also observed @. Me- 
chanical effects of light in optical resonators were studied 
in while cavity- assisted nondestructive laser cooling of 
atomic qubits was analyzed in [7] . A laser cooling method 
that can be used at large detuning and low saturation to 
cool particles inside an optical cavity was proposed in 
Ref. @- A significant speed-up of the cooling process 
was found in [9| while fast cooling of trapped ions using 
the dynamical Stark shift was described in (Toj . 

Via engineering superconducting elements as artificial 
atoms and coupling them to a photon field of a resonator 
or to vibrational states of a nanomechanical resonator 
one can demonstrate interesting related phenomena such 
as single artificial atom lasing or cooling. In particu- 
lar, schemes to ground-state cooling of mechanical res- 
onators were proposed in fllj j . A flux qubit was experi- 
mentally cooled [12j | using techniques somewhat related 
to the well-known optical sideband cooling methods (see, 
e.g., Ref. [l[ and references therein). Lasing effects of a 
Josephson-j unction charge qubit, embedded in a super- 
conducting resonator, was experimentally demonstrated 
in [13j . Single-qubit lasing and cooling at the Rabi fre- 



quency was proposed in [Tij , while a mechanism of simul- 
taneously cooling of an artificial atom and its neighboring 
quantum system was analyzed in [15|. Few-qubit lasing 
in circuit QED was discussed in Ref. [16[ . A LC oscillator 
can be cooled via its nonlinear coupling to a Josephson 
flux qubit [17] . The cooling of a nanomechanical res- 
onator via a Cooper pair box qubit has been recently 
suggested in Ref. [l8| while cooling carbon nanotubes to 
the phononic ground state with a constant electron cur- 
rent was achieved in [L9| . Further interesting works on 
cooling micro- and nanomechanical resonators were pre- 
sented in Refs. US-Hi]. 

Here, we describe a cooling scheme via coupling a 
pumped multiparticle ensemble (i.e. artificial atoms or 
qubits) to a single mode of a quantum LC circuit (see 
Fig- Hi- Our motivation is to present an efficient method 
allowing for a rapid cooling of the resonator mode. The 
multiqubit system can be formed by an independent N- 
particle sample or by collectively interacting N parti- 
cles. By independent, we mean that each particle spon- 
taneously decays individually and all of them are maxi- 
mally coupled with the oscillator mode and with the same 
phase. Collectively interacting particles means that their 
interactions are mediated by the environmental electro- 
magnetic field reservoir such that their decay is of a col- 
lective nature. In this case, the particles are close to each 
other on a scale smaller than the emission wavelength and 
coupled with the same strength to the quantum oscillator 
mode. The advantages or disadvantages regarding the in- 
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FIG. 1: An ensemble of TV independent three-junction flux 
qubits is coupled to a LC quantum circuit by their mutual 
inductances. An AC magnetic flux <3?(i) drives the qubits. 
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terparticle interactions to the cooling phenomena of the 
quantum oscillator degree of freedoms will be discussed 
in detail. In particular, we found that the cooling phe- 
nomenon is better for independently interacting qubits if 
the quantum dynamics of the LC oscillator is slower than 
that of the qubits. However, the cooling effects may oc- 
cur faster for collectively interacting qubits. Apart from a 
fundamental interest, these systems have a great feature 
in various applications such as novel quantum sources 
of light (single photon sources, for instance), quantum 
processing of information or entanglement. However, at 
MHz frequency ranges thermal fluctuations affect con- 
siderably the LC oscillators, i.e., populate their energy 
levels and induce additional dccohcrences. Therefore, a 
suitable method to cool these systems can be very useful. 

The paper is organized as follows. In Section II, we 
introduce the system of interest. Section III describes the 
obtained results. We finalize the article with conclusions 
presented in Section IV. 

II. APPROACH 

We describe the cooling effects of a quantum oscillator 
mode, i.e. a quantum LC circuit coupled inductively with 
a collection of two- level Josephson flux qubits (see Fig. [I]). 
The two-level particles are pumped with a moderately in- 
tense magnetic flux and damped spontaneously via their 
interactions with the environmental electromagnetic field 
reservoir. The single-particle spontaneous decay rate is 
7. Both subsystems interact with thermostats at effec- 
tive temperatures T\ and T2. We shall consider that the 
particles are independent or collectively interacting. The 
frequency of the oscillator is much lower than the qubit's 
tunnel splitting, i.e. lj c <C A. Therefore, the qubit is 
driven with Rabi frequencies near resonance with the os- 
cillator frequency that affect the oscillator, increasing its 
oscillation amplitude |30J. Near the symmetry point (i.e. 
the energy bias e between the flux states is small) and 
after transformation to the qubit's eigenbasis, the Hamil- 
tonian describing the multiqubit systems is: 

N 

H = AEa Z i/2 + flp sin29a X i cos (wt)} + gj c a^a 

i=\ 
N 

— g (cos 2#er Z i — sm29a X i)(a + aJ), (1) 
»=i 

where the first term describes the qubits, each with 
the transition frequency AE = \/A 2 + e 2 , while the 
second one considers their driving by an applied AC 
magnetic flux with amplitude f2o = f2o sin 29 and fre- 
quency uj. Here, cot 29 = e/A with cos 2$ = e/AE and 
sin 29 = A/AE. The last two terms describe the oscil- 
lator with frequency lu c = 1/yLC as well as the qubit- 
oscillator interaction, respectively. Here, g ps MI p I c0 , 
where M is the mutual inductance, I p the magnitude of 
the persistent current in the qubit, and I c q = \J lo c /2L 



the amplitude of the vacuum fluctuations of the current 
in the LC oscillator, and a are the creation and anni- 
hilation operators corresponding to the oscillator degrees 
of freedom, while a a (a € {x,y,z}) are the Pauli ma- 
trices operating in the dressed flux basis of the qubit 
subsystem. As A » w c , the transverse coupling in the 
Hamiltonian (fTJ) is transformed into a second-order lon- 
gitudinal coupling by employing a Schrieffer- Wolff type 
transformation [lj, [3l| , i.e. Us = exp (iS) with 

N 

S = (g/AE) sin26»(a + a f ) ^ o yi . 

i=i 

By further using the rotating wave approximation with 
respect to uj and diagonalizing the qubit term as well as 
applying the secular approximation, i.e. omitting terms 
oscillating with the generalized Rabi frequency, one ar- 
rives at the following Hamiltonian describing the interac- 
tion between the multiqubit system and the LC oscillator: 

JV 

H = iJ c a^a + ^{nR zi /2 + g(R { l ) _a + a^R ( i ) + ) 

i=l 

+ g Q {aa) +a)a)R zl /2}. (2) 

Here, we have further assumed that the generalized Rabi 
frequency is of the order of Li c , that is ps cj c . In 
Eq. ©, g=g cos 29 sin 2£ gives the qubit-oscillator cou- 
pling strength while <?o=2g 2 cos 2£ sin 2 29 j AE accounts 
for a small frequency shift of the qubit's frequency. Fur- 
ther, 

cot 2^ = 5lj/£Io, (3a) 
cos 2 C = [1 + 5u/Q}/2, (3b) 
sin 2 £ = [1 - 5u/fi\/2. (3c) 

where 5oj = AE — uj and where 17 = \J (5lo) 2 + O 2 , 
stands for the generalized Rabi frequency. The dressed- 
state qubit operators R^l = \a)a{j3\ describe the in- 
ternal transition in the ith particle between the dressed 
state |/3) and \a) for a ^ /3 and population for a = /3, 
{a, ft G +, — }, and obey the standard commutation rela- 
tions of su(2) algebra, that is 

[R$»R%] = h{^>R% - s p , a R%), 

where a,/? € { + , — }. R Z i = — is the dressed- 
state inversion operator for the ith particle. 

In the mean-field, dipole, Born-Markov and secular ap- 
proximations, the combined system is characterized by 
the following master equation: 

—p + i[H,p] = -AaP- A c p. (4) 

The quantum dissipation due to spontaneous emission 
into surrounding electromagnetic field reservoir is de- 
scribed by the A a p term which for N independent qubits 
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can be represented as follows: 



JV 



7 .[R%Rf_p]} + H.c. 



(5) 



For N nonindependent radiators, i.e. for collectively in- 
teracting particles, the corresponding damping is: 



/v 



AaP = fro ^ ' R *3P\ + ^+ l R +- < R ~+P] 

+ 7 _[i? ( ! ) + ,i?i J lp]}+i?.c. (6) 
The damping rates are given by the following expressions: 



7+ = 7 cos 4 £, 

7- = 7 sin4 & 
70 = 7 sin 2 2£/4. 



(7a) 
(7b) 
(7c) 

The last term in Eq. Q characterizes the damping of the 
quantum oscillator mode and is given as follows: 

A c p = k(1 + n(w c )) [a', ap] + ™(u c )[a, a'p] + H.c. (8) 

Here, n(oj c ) is the mean thermal photon number corre- 
sponding to the resonator frequency tu c while k is the res- 
onator decay rate. We have omitted the coherent part of 
the dipole-dipole interaction in Eq. which is justified 
if the Rabi frequency dominates over the dipole-dipole 
induced energy shifts. 

A general analytical solution of Eq. Q is not evident. 
However, one can obtain its solution for different regimes 
of interest, namely in the bad or good cavity limit. There- 
fore, in the next Section, we proceed by investigating the 
properties of Eq. (0J when the qubit's quantum dynamics 
is faster than the one of the quantum oscillator, i.e. in 
the good cavity limit (321 ] . 



III. RESULTS AND DISCUSSIONS 

We assume a moderately intense pumping field, i.e., 
^ ^ {7jSV^V} and a high quality resonator such that 
7 » gV~N ^> n for N independent particles or Q ^> 
{Nj,gy/N}, and N 7 > gy/N > k for N collectively 
interacting particles. Therefore, in this case, the qubit 
subsystem achieves its steady-state on a time scale faster 
than the resonator field and, thus, the qubit variables 
can be eliminated to arrive at a master equation for the 
resonator field mode alone: 



P : 



T-{a'ap — apa'} — T + {aa* p — a' pa} + H.c, (9) 



where an overdot means differentiation with respect to 
time and 



r_ = n(l + n{u>cj) 
r+ = Kn(u c ) + A. 



B, 



(10a) 
(10b) 



The physical meaning of the parameters in Eq. ([9]) is as 
follows: r + (r_) describes the process of increasing (de- 
creasing) of the photon number in the resonator mode. 
The interplay between r + and T_ leads to lasing or cool- 
ing of the quantum LC circuit. 

For an independent N particle system one has: 
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(#4 



and B 



r_L 



(R- 



where 



r_L = 470 + 7+ + 7- 



(11) 



(12) 



We have considered here that (i?^!) are identical for 



*/3/ 



all i £ {1,2,- ■■ ,N} and, thus, (R^) = (R a0 ) where 
{a, P € +, — }. The expectation values for {R a p) are cal- 
culated in the absence of the resonator mode. Therefore, 
from Eq. ([5]), we have: 



(R, 



7- 



7- +7^ 



and (R-) 



7+ 



7- + 7-f 



(13) 



The expressions (fTTj) are valid for any N satisfying the 
restrictions imposed in the beginning of the Section in- 
cluding N = 1. 

For an ensemble of collectively interacting N particles, 
we obtain the following relations for A and B: 

~2 ~2 

A = ^-(R + _R_ + ), and B = £-(R-+R+-), (14) 

1 _L 1 _L 

where the collective decay rate is given as follows: 

f±=rx + (7--7+)(^}- (15) 

One can observe here that the decay rate T± has a 
contribution arising from all particles, i.e., the term 
proportional to (R z ). Here, in contrast to indepen- 
dent qubits, collective operators were introduced, that 
is, R a p = J2iLiRal- Note that to obtain Eqs. (fTTl) . 



we decoupled the involved multiparticle correlators- an 
approximation valid for larger N, i.e., JV > 1. How- 
ever, the corresponding expressions for N = 1 are iden- 
tical to Eqs. (fTTj) but with the single-particle decay rate 
instead of collective ones. The steady-state expecta- 
tion values for the collective correlators entering into the 
above expressions, i.e. Eqs. (|14p. can be estimated from 
the steady-state solution of the master equation Eq. ([6|) 
describing the strongly driven particles in the absence of 
the resonator [33, kBJI : 



(Rz) = 



N(l-f 2 + N ) + f(N + 2)(f N -1) 



(R+-R-+) 
(R-+R+-) 



1-/ 
/ 

1-/ 



(Rz), 



(16a) 
(16b) 
(16c) 
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FIG. 2: (color online) The mean photon number (n) into the 
quantum circuit as function of 8ui and different numbers of 
independent qubits. The solid line is for TV = 1, the long- 
dashed line stands for TV = 10, while the short-dashed curve 
corresponds to TV = 30. The dotted curve shows the satura- 
tion photon number no for TV = 30 qubits. Here, n(u c ) — 4, 
A/2tt = 3- 10 9 Hz, e = 0.01A, g/2-rv = 10 6 Hz, uj c /2tt = 10 7 Hz, 
7 /2tt = 10 5 Hz, k/2tt = 10 3 Hz and Q = \/^ 2 ~(5lo) 2 . 



FIG. 3: (color online) The mean photon number (a* a) of the 
quantum oscillator as the function of Slj and different numbers 
of collectively interacting qubits. The solid line is for TV = 10 
while the long-dashed line stands for TV = 15. The short- 
dashed curve corresponds to TV = 30 while the dotted one to 
TV = 30 but k/2tt = 10 2 Hz. Other parameters are the same 
as in Fig. ©. 



where / = 7+/7-. For TV = 1 one obtains Eqs. (p~3j) . 

The Eq. ([9]) has an exact steady-state solution. For 
instance, the steady-state expectation values for the di- 
agonal elements of Eq. ^ can be obtained from the re- 
lation: 



exp| 



-aa* a] 



(17) 



where a = m(7y) with i] = r_/r+ and Z is determined 
by the requirement Tr{p} = 1. Evidently, the expec- 
tation values of the operators needed for evaluating the 
properties of the quantum oscillator are obtained from 
Eq. (fTTj) . In particular, the oscillator mean photon num- 
ber, i.e. (n) = (a^a) = Tr{a^ap}, and its second-order 
correlations can be determined from the following expres- 
sions: 



(a f a) = 



1 



ry-1' 



a 2 ) = 



fa l) 2 



(18a) 
(18b) 



to such an extent that the photon second-order correla- 
tion function, i.e. 

J2W _ ( afa2 ) 



ff w (0) 



(oto) 



equals with 2 which means that the photon statistics is 
always super-Poissonian. Note that the mean photon 
numbers obtained with the help of Eq. © or Eq. (fTT)) 
should be below the photon saturation number uq which 
for TV interacting qubits reads approximately as: 



n o =r ± ( 7 ++7-)/(ff 2 A0 



(19) 



Fig. ([2]) depicts the mean photon number in the oscil- 
lator mode which is coupled with TV independent qubits. 
We have used typical parameters here (see, for instance, 
30]). To elucidate the role of many particles regard- 
ing the cooling issue, we fix the involved parameters and 
change the number of qubits. Already for TV = 10 par- 
ticles, the cooling efficiency is significantly improved in 
comparison to the single-qubit case, i.e. TV = 1. Better 
cooling can be achieved, that is (n) <C n(ui c ), by increas- 
ing further the number of qubits (see the short-dashed 
curve in Fig. [3]). Evidently, the qubits are in their lower 

dressed-state when cooling occurs, i.e. (R ) > (R ++ ). 

The diagram showing the energy levels of the qubit and 
oscillator indicating the cooling cycle with photon emis- 
sion/absorption can be found in Refs. [25l |26||. 

Further, we turn to cooling effects via collectively inter- 
acting particles. Fig. ([3]) shows the mean photon number 
in the quantum oscillator mode as the function of 5u>. 
The mean photon number (n) is well below the ther- 
mal mean photon number fi(ui c ), however the cooling is 
not so significant as for independent qubits (compare the 
short-dashed curves in Fig. [2] and Fig. [3]). The reason 
is that the decay rate T± is dependent on the number 
of qubits and, thus, (n) = r + /(r_ — T + ) is smaller 
than the corresponding one for independent qubits since 
T_ -r+ = n-g 2 (R z )/f ± where (R z ) is given by Eq. (JTSJ) 
(in other words, for collectively interacting particles, we 
do not have a factor TV in the denominator). However, 
adjusting the involved parameters, one can improve the 
cooling efficiency in general (see the dotted line in Fie 
Note that the coupling of qubits can be controlled [3£ 

Finally, we discuss the time scaling for the cooling phe- 
nomenon. We observe that cooling rates depend on the 
number of qubits and, therefore, the cooling may occur 
faster in both schemes. However, the faster decay rate in 
our approach is the qubit spontaneous emission. Thus, 
the cooling phenomena can not occur faster than 7 _1 for 
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independent qubits or (A^) -1 for collectively interacting 
particles, respectively. Therefore, in general, the cooling 
processes are faster for collectively interacting particles. 

IV. SUMMARY 

In summary, we described a scheme that is able to cool 
a quantum LC circuit coupled inductively to externally 
pumped artificial particles (Josephson flux qubits) and 



damped through their interaction with the environmen- 
tal electromagnetic field reservoir. The qubits may in- 
teract collectively or they are independent. If the qubits 
dynamics is faster than the one of the LC oscillator, the 
cooling of the oscillators degrees of freedom occurs when 
controlling the qubit quantum dynamics. We found that 
the cooling phenomenon is better for an ensemble of in- 
dependent qubits. However, in general, the cooling pro- 
cesses are faster for collectively interacting particles. 
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